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We investigate the dimensionality of inhomogeneous spin modulation patterns in the cuprate fam- 
ily of high-temperature superconductors with particular focus on 1/8-doped lanthanum cuprates. 
We compare one-dimensional stripe modulation pattern with two-dimensional chackerboard of spin 
vortices in the context of nuclear quadrupole resonance(NQR) and muon spin rotation(/iSR) ex- 
periments. In addition, we also consider the third pattern, a two-dimensional superposition of spin 
spirals. Overall, we have found that none of the above patterns leads to a consistent interpreta- 
tion of the two types of experiments considered. This, in particular, implies that the spin vortex 
checkerboard cannot be ruled out on the basis of available NQR//iSR experimental results. 

PACS numbers: 



I. INTRODUCTION 

The apparent connection between superconductivity 
and magnetism in cuprates, pnictides and other uncon- 
ventional superconductors is not, at present, well under- 
stood theoretically. A possible mechanism behind this 
connection is through the formation of slowly fluctu- 
ating inhomogeneous magnetic structures emerging as 
a result of strong tendency towards electronic phase 
separation [lj. Electronic phase separation fluctuating in 
space and time is difficult to investigate experimentally. 
In this respect, the lanthanum family of high temper- 
ature cuprate superconductors at the doping level 1/8 
has long been the subject of intense interest, due to 
the fact that electronic spin and charge modulations in 
these materials are static in time and periodic in space. 
The character of these modulations, however, is still the 
subject to mutually exclusive propositions. The pur- 
pose of the present paper is to examine these propo- 
sitions from the viewpoint of nuclear quadrupole res- 
onance (NQR) experiments and muon spin rotation 
(/iSR) experiments 

At present, the dominant viewpoint is that the 
spin and charge modulations in 1/8-doped lan- 
thanum cuprates have the character of one-dimensional 
stripes Motivated by the STM observations of 

checkerboards one of us previously investigated the 

two-dimensional alternatives to the stripe interpretations 
— initially in the form of grid (l3l - [T5| and then, after neu- 
tron scattering evidence against grid appeared (Itf, in the 
form of spin vortex checkerboard (l7l TLsj . Refs.[l3l li~5l . 
[l7l [Lsj also contain detailed discussion of the arguments 
in favor and against the stripe proposition. The read- 
ers are referred to Refs. [l9i - [27| for various considerations 
pertaining to the formation of spin vortices and other 
non-collinear spin superstructures. 

Static spin modulations create the distribution of local 
magnetic fields, whose effect is observable by both NQR 
and /iSR experiments. The advantage of these experi- 
ments is that one can make experimental predictions for 
a given spin modulation pattern essentially without fur- 



ther theoretical assumptions. The disadvantage is that 
these techniques do not resolve the wave vectors of spin 
modulations. 



The comparative analysis of the stripe and grid super- 
structures from the viewpoint of NQR experimcntswas 
done in Refs.fl5| and from the viewpoint of /iSR in 
Ref. Q . The results of the NQR analysis indicated that 
the stripe pattern would lead to significant discrepancies 
with experiment due to a certain singularity in the dis- 
tribution of local magnetic fields. At the same time, the 
grid pattern did not lead to that singularity and, as a 
result, was better compatible with NQR experiments. In 
general, however, the NQR experiments were even better 
compaiblc with a totally random Gaussian distribution 
of the local fields, which, in turn, would be incompatible 
with the static modulated spin response observed by elas- 
tic neutron scattering 0,0. The mSR results @, on the 
other hand, were more consistent with the stripe pattern 
than with the grid precisely because of the presence of 
the above mentioned singularity on the distribution of lo- 
cal magnetic fields. At the same time, it was not possible 
to reproduce the /iSR results quantitatively from the full 
three-dimensional stacking of either stripes or grid. 

The above rather incoherent phenomenology is further 
complicated by the discrepancy in the estimate of the 
spin modulation amplitude from /iSR (about 0.3/is) and 
from neutron scattering (0.1 — 0.15/is). The NQR exper- 
iments are more consistent with the latter estimate. 



So far, the attempts to analyze NQR/^SR data have 
not included the spin vortex lattice. It is our goal in this 
paper to check whether this pattern may help to reconcile 
NQR/ yuSR/neutron scattering phenomenology. For com- 
pleteness of the analysis, we have also done NQR/^SR 
calculations for a coherent superposition of two orthog- 
onal spiral harmonics previously mentioned in Ref. [171]. 
We further investigated possible muon sites which are 
different from the "standard" site identified in Ref. 1281 . 
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FIG. 1: Site-centered, commensurate spin superstructures (a) stripes, (b) spin-vortex lattice, and (c) 2D superposition of spin 
spirals. 



II. THEORETICAL FORMULATION 

A. Spin Structures 

In 1/8-dopcd lanthanum cuprates, the primary ex- 
perimental evidence of spin modulations is based 
on the four-fold splitting of the antiferromagnetic 
± 6=Z-, ~ ± S^f) j were a the in-plain lattice period 
and 5 ps 1/8. This splitting indicates that the antifer- 
romagnetic order is modulated with period a/5. The 
experiments also indicate that spin polarizations mainly 
lie "in-plain" . All possible modulations that may be con- 
sistent with the above four peaks can be represented as 



S'i x cos(qi • Tij + ip lx ) 
5i y cos(qi • + ip ly ) 

S 2x cos(q 2 ■ + ip 2x ) 
S 2y cos (q 2 • + ip 2y ) 



, (1) 



where i and j are the indices of the square lattice, 
Tij are the positions of the lattice sites, S^- the static 
spin polarization values for each site, qi = (8— , 0) 
and q2 = (0, 5—) are the modulation wave vectors, 
{Six,Si v ,S 2x S 2y } and {ipi x , ip ly , (p 2x , ¥>2y} are respec- 
tively the amplitudes and the phases of four possible lin- 
ear polarization harmonics. Below, we use variable So 
to denote the maximum spin polarization of a given pat- 
tern. Wc also use Bohr magneton units for electronic 
spin polarizations, so that the maximum polarization of 
spin 1/2 corresponds to S = Hb- 

Possible specific patterns of general form([T])include: (i) 
stripes, e.g. {S lx = 0,5^ = S ,S 2x = 0, S 2a = 0} 
[Fig. ma)]; ( n ) spirals, e.g. {S lx = S ,S ly = S ,S 2x = 
0,S 2y = 0,(pi x - f\y = 7r/2}; (iii) grid, e.g. {S\ x = 
0, Siy = S /2,S 2x = 0,S 2y = S /2}; (iv) spin vor- 
tex checkerboard, e.g. {Si x = 0,Si y = So/V2,S 2x = 
So/V2, Siy = 0} [Fig.QJb)]; (v) Superposition of two or- 
thogonal spirals, e.g. {S\ x = Sq/2,Si v = S /2,S 2x = 
S / y/2, S 2y = S / v2, (fix - fiy = tt/2, ip 2x - ip 2y = tt/2} 
[Fig. [He)]. All these patterns are locally stable in the 
approximation of staggered spin polarizations [T^. The 
stripe and the spiral interpretations are one-dimcsional 



and, therefore, require two kinds of domains with orthog- 
onal modulations to account for the four-fold splitting of 
the antiferromagnetic peak. 

Neutron scattering study of Ref. [l6| indicates that the 
above spin modulation harmonics should be transversely 
polarized, which implies that S\ x = and S 2y = 0. 
This, in turn, is inconsistent with the spiral patterns (iii) 
and (v) and with the grid pattern (iii), thus leaving the 
transversely polarized stripes (i) and spin vortex super- 
lattice (iv) as the only two remaining propositions. The 
comparative analysis of these two propositions is the pri- 
mary goal of this work. Below, however, we also present 
NQR//iSR calculations for the two-spiral pattern (v), 
which were done, in part, to make the analysis complete 
and, in another part, because as obvious from Fig. [He)] 
this pattern appears to be consistent with the minimum 
of the pseudogap along the diagonal lattice direction [2^]. 

When necessary for the /xSR calculations, the three- 
dimensional arrangement of the above spin modulations 
was chosen as follows. For stripes, we use the orthogo- 
nal stripe directions in the adjacent layers, and in sec- 
ond adjacent layer — half-a-period parallel shift. Spin 
vortices checkerboards and the two-spiral superstructures 
are both shifted in the adjacent layers by approximately 
two lattice periods in both the a— and 6-directions, so 
that the spin-poor regions in one layer are located on 
the top of the spin-rich regions in the adjacent layer. 
In the case of stripes, the above arrangement is what is 
proposed by Tranquada, even though theoretical doubts 
about such proposition exist [HI, [l8|. For the two other 
superstractures, the above arrangement is based solely on 
the basis of Coulomb repulsion between the same-charge 
regions, and the assumption that spin-poor regions accu- 
mulate charge carriers and thus positively charged, while 
spin-rich regions are negatively charged. 

In the following, we also explore the noisy modifica- 
tions of the above modulations. The spin noise is to be 
uncorrelated between different lattice sites. For each site, 
it will have two in-plane components, each randomly se- 
lected from the Gaussian distribution with zero average 
and root-mean-squared deviation denoted as AS. The 
noise was incorporated in the NQR calculations by an- 



3 



alytically modifying the formulas presented in the next 
section, while the /j,SR calculations were based on the 
actual sampling of the three dimensional spin structure 
with or without the noise. 



B. NQR lineshapes 

At low temperatures, static spin modulations lead to 
a broad distribution of local hyperfine fields that affect 
NQR frequencies and transition rates. In this paper, we 
attempt to describe the experimental low-temperature 
63 Cu NQR lineshape reported in Ref.[||. Nuclear iso- 
topes 63 Cu have spin 3/2. The NQR Hamiltonian for 
each nuclear spin is assumed to have form 

H = ^|4 2 -i/(/ + l)|- 7 Cu^H.I (2) 

where I is the nuclear spin operator, h is the Planck's 
constant, 7c u is the gyromagnctic ratio, H is the local 
hyperfine field created by static spin modulations, vq is 
the quadrupolar parameter, and the z-axis is perpendic- 
ular to the CuC>2 plane. The hyperfine local field at the 
positon Vij is given by 

NN 

H(r y ) =AS ij + B Y, S kh (3) 

kl 

I 



where A = 38 [kOc/jUe] and B — 42 [kOc/^ B ], and the 
sum is taken over the four nearest neighbors of spin S^-. 

The calculation of the NQR lineshape (echo intensity 
as a function of the frequency of NQR echo pulses) re- 
quires computing six transition frequencies between the 
four eigenstates of the Hamiltonian ([2]), weighing these 
transitions by the square of the appropriate matrix ele- 
ment and by the square of the frequency itself and then 
averaging over all possible values of the hyperfine field 
associated with a given spin superstructure. The frame- 
work and the NQR parameters of our calculations includ- 
ing the A- and B- NQR lines and broadening of these 
lines are identical to those used in Refs.jH, El, where 
stripe and grid superstructures were considered. Here 
we only describe the last step associated with comput- 
ing the hyperfine field distribution due to the spin vor- 
tex checkerboard [Fig. GXb)] and the two-spiral modula- 
tion[Fig.QJc)]. We assume that the periods of both mod- 
ulations are very close but not exactly equal to 8 lattice 
periods. As a result, both modulations are incommen- 
surate, and the resulting distribution of hyperfine fields 
involves the averaging over the "ergodic" phase of these 
modulations. 

For the spin vortex checkerboard, the resulting distri- 
bution of hyperfine fields can than be obtained analyti- 
cally: 



p{H) = -2- x < 



d{ 



(C 2 -i)C 



C 2 +i— 



(H/Hr, 



(C2 - 1) ( C2 -€) ( C2 + 1 -€) 



ifO<F<^ 



if % < H < H, n 



V2 



elsewhere. 



(4) 



Here Nq is the normalization constant, and H m = 
S [B(2+V2)-A) is the maximum hyperfine field due in- 
duced by the spin vortex checkerboard. This distribution 
has Van Hove singularity at H = H m / a/2 ■ It is shown in 
Fig. [2] (a, left column, blue line) and compared to the dis- 
tribution for the stripe structure (magenta line), which 
also has a Van Hove singularity but at H = H m . 

The distribution of the local hyperfine fields for the 
two-spiral superstructure is identical to that of the stripe 
superstructure. 



C. /iSR signal 

The time evolution of the muon polarization vector P M 
for polycrystalline samples in the approximation of static 
local fields and without additional depolarization effects 



can be calculated as follows 



1 2 f°° 
P » W = 3 + 3 J Q P{H) C ° S t] dA 



(5) 



where p(H) is the probability distribution for the ab- 
solute value of the local field Hat the position of the 
muon and 7 M = 851.61 [MHz/T] is the muon gyromag- 
netic ratio. We units of are to be normalized such that 
-F^(O) = 1. We perform the calculations for the muon site 
identified in Ref. HU to be = (0.253,0,0.162) [in the 
relative lattice units]. 

We obtain the distribution p{H) by sampling dipolar 
fields 



H 



E 



3(S« 



..2 ( 



(6) 



where r,,- is the displacement vector between the muon 
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FIG. 2: Left column: Distribution of hyperfine fields for an incommensurate spin lattice. Middle column: NQR simulation 
compared to the experimental data of Hunt et al. 0. Right column: /^SR simulation compared to the experimental data of 
Luke et al. (a) So = 0.32, AS = 0, (b) So = 0.20, AS = 0, and (c) So = 0.20, AS = 0.5 S . For the left and middle column 
the curves for the 2D spin spirals coincide with the curves for spin stripes. 



site and the copper positions on the lattice. The sampling 
was done for all possible muon sites in a cluster of 32 x 
32 x 16 spins. 



III. RESULTS AND DISCUSSION 
A. Joint analysis of NQR and /iSR experiments 

The results of our calculations of NQR lincshapes and 
/iSR signals are presented in Fig. [2] The first column of 
plots contains the distribution of hyperfine fields for the 
calculation of the NQR lineshape. The second column 
represents the NQR lineshapes themselves. The third 
column contains the corresponding [iSR signals. 

The first row of panels in Fig. [5] labeled as (a) rep- 
resents calculations for the modulation amplitude Sq = 
0.32/is and no spin noise, i.e. AS = 0. This choice is 
most appropriate for fitting the initial behavior of /iSR 
signals even though significant discrepancies with exper- 
iment remain at intermediate and long times. The NQR 



lincshapes in Fig. HJa) for all modulation patterns con- 
sidered exhibit apparent disagreement with experiment. 
The experimental lineshape contains a single peak, which 
is likely associated with the van Hove singularity in the 
function converting the hyperfine fields into the NQR fre- 
quencies. This van Hove singularity is to be contrasted 
with the the multiple theoretical peaks appearing in the 
middle panel of Fig. HJa) and originating from the van 
Hove singularities in the distribution of the hyperfine 
fields themselves [see the left panel of Fig. dja)]. 

The second row of panels in Fig. [2] labeled as (b) rep- 
resents our best attempt to fit the NQR lineshape by 
reducing the amplitude of spin modulations, still with- 
out introducing spin noise. The optimal modulation am- 
plitude was Sq = 0.2/is- For this value, the principal 
peaks originating from the van Hove singularities shown 
in the left panel of Fig. [U[b) coincide with the peak of of 
the function converting the hyperfine fields into the NQR 
frequencies. The experimental NQR line is still notice- 
ably broader then cither of the theoretical predictions. 
The corresponding theoretical /iSR signals arc noticeably 
slower than the experimental one. 
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Finally, the third row of panels in Fig. [2] labeled as (c) 
represents our best attempt to further improve the agree- 
ment with the experimental NQR lincshape by adding 
spin noise: So = 0.2/is, AS = 0.5So- In this case, all 
three patterns considered lead to a good agreement with 
the experimental NQR shape, but a significant discrep- 
ancy with the experimental /iSR results remains. 

The overall outcome of our calculations is that nei- 
ther of the patterns considered can accurately account 
for both the NQR and the /J.SK experiments with the 
same values of So and AS. Either of the patterns can 
well reproduce the NQR lineshape with a smaller modu- 
lation amplitude and in the presence of a strong degree 
of spin noise. The /iSR experiments are best fitted with 
a larger modulation amplitude and in the absence of the 
spin noise. Still, even the best fits of the fiSR experiments 
exhibit significant disagreements with the experiments. 

Our experience suggests that the origin of the above 
persistent disagreement in the shape of the /xSR sig- 
nal is beyond the specific choice of the spin modula- 
tion pattern. We expect a significant disagreement for 
any pattern of the static local spin configurations. It 
originates from the fact that the three-dimensional ar- 
rangement of static spin polarizations unavoidably leads 
a three-dimensional distribution of magnetic dipole fields 
experienced by muons. For a three-dimensional distribu- 
tion, the probability of the zero absolute value of mag- 
netic field is normally equal to zero. This probability 
is proportional to the Fourier transform of Py.(t) — 1/3 
at zero frequency. Therefore, zero probability for the 
zero absolute value of the local magnetic field implies 
that J Q [Pf_i(t) — 1/3] dt = 0. The experimental /iSR sig- 
nal used for comparison with theory in Fig. [5] apparently 
does not satisfy the above condition. It is not clear to us 
whether the above discrepancy has something to do with 
the inadequacy of the theoretical assumption of static lo- 
cal fields, or with experimental uncertainty in terms of 
identifying the 1/3-lcvcl of the /xSR signal. 

Irrespectively of the above disagreement, it was argued 
by Kojima et al.Q, that large amplitude of /iSR oscilla- 
tions indicate that the spin structure is one-dimensionally 



modulated. This conclusion was based on the the com- 
parison between the results for one-dimensional stripes 
and two-dimensional grid. However, our calculations 
indicate that, in the case of the two-dimensional spin- 
vortex checkerboard, the resulting oscillations are even 
larger than for the stripe pattern. This indicates that 
large amplitude of fiSK oscillations as such cannot dis- 
criminate between stripes and spin vortices and hence 
between ID and 2D spin patterns. 

Finally, we have attempted to explore the possibility 
that the muon site is different from the one identified 
by Hitti et al.[28j. Trying the muon positions in the 
ac-plane lA away from the apical oxygen, we were able 
to improve the /xSR fit for the stripes and 2D spirals 
with muon position at = (0.110,0.0,0.114) and with 
So = 0.15/iBj AS* = 0. In this case, the NQR spectrum 
has a single main peak, which is much narrower than the 
one observed experimentally, but with the maximum at 
the experimentally observed frequency. No comparable 
improvement of the /iSR fit was found for the spin vortex 
checkerboard. 



IV. CONCLUSIONS 

The outcome of the analysis presented in this paper 
indicates that spin vortex lattice leads to the singular- 
ity in the distribution of local magnetic fields, which is 
quite similar to the one produced by the stripe pattern. 
As a result no clear qualitative differences emerge be- 
tween the NQR//xSR predictions based on stripes or spin 
vortices, while, in the both cases noticeable quantitative 
discrepancies with the experimental results remain. Sim- 
ilar discrepancies also appear for the 2D superposition of 
spin spirals. It is only clear that the observation of the 
large amplitude of fiSK oscillations does not amount to a 
qualitative argument to support stripes against the spin 
vortex checkerboard. 
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